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We present first principles calculations of the magnetic moments and magnetic anisotropy energies 
of small Fe, Co and Ni clusters on top of a Ag(100) surface as well as the exchange-coupling en- 
ergy between two single adatoms of Fe or Co on Ag(100). The calculations are performed fully 
relativistically using the embedding technique within the Korringa-Kohn-Rostoker method. The 
magnetic anisotropy and the exchange-coupling energies are calculated by means of the force the- 
orem. In the case of adatoms and dimers of iron and cobalt we obtain enhanced spin moments 
and, especially, unusually large orbital moments, while for nickel our calculations predict a complete 
absence of magnetism. For larger clusters, the magnitudes of the local moments of the atoms in the 
center of the cluster are very close to those calculated for the corresponding monolayers. Similar to 
the orbital moments, the contributions of the individual atoms to the magnetic anisotropy energy 
strongly depend on the position, hence, on the local environment of a particular atom within a 
given cluster. We find strong ferromagnetic coupling between two neighboring Fe or Co atoms and 
a rapid, oscillatory decay of the exchange-coupling energy with increasing distance between these 
two adatoms. 

PACS numbers: 71.24.+q, 71. 70. Cm, 75.30.Gw, 75.30.Hx 



I. INTRODUCTION 

Magnetic nanostructures such as impurities, clusters 
and wires on top or in the uppermost layers of surfaces 
are of special interest for nano-scale technologies, in par- 
ticular, regarding their possible application as magnetic 
nano-devices and high-density magnetic recording me- 
dia. A quantitatively correct description of the mag- 
netic properties of such structures, namely, the magni- 
tude and the orientation of spin and orbital moments, 
magnetic anisotropy energies and the magnetic interac- 
tions, is, therefore, an important issue to be addressed. 
Concomitantly, the understanding of the changes of phys- 
ical properties from nanostructures to thin films or even 
bulk systems has always been a fascinating theoretical 
challenge. 

Because of the lack of translational symmetry tight- 
binding (TB) methods have been an efficient tool to 
study larger clusters. By using a tight-binding Hub- 
bard Hamiltonian in the unrestricted Hartree-Fock ap- 
proximation, Pastor and co-workers revealed the size and 
structural dependence of magnetic properties of free Cr n , 
Fe„ and Ni„ (n < 15) clustersa, and also the exchange in.- 
teraction and local environments effects in Fe n clusterso. 
By including a spin-orbit coupling term into the Hamil- 
tonian, they also investigated various effects on the mag- 
netic anisotropy energy (MAE) of small unsupported Fe 
clusterstJ and, recently, of Co n clusters on Pd(lll)H Fi- 
nite temperature magnetism of small clusters, remark- 
ably different from that of bulk systems, has also been 
studied in terms of a similar approach by takingr.into 
account both electronic and structural excitationsQ. A 
great advantage of the TB methods seems to be that 



they easily can be combined with molecular dynamics 
calculations enabling thus investigations of relaxation ef- 
fects which provecLto be important—in determining the 
magnetic momentsTB and the MAEo of transition metal 
clusters. 

The embedding technique based on the Korringa- 
Kohn-Rostoker (KKR) Green's function method in the 
local spin-density approximation (LSDA) has been ap- 
plied to the magnetism of trapsitipn metal adatoms and 
clusters deposited on surfacesBB The main feature of 
this approach is that the interaction between adatoms 
and host surface atoms can be analyzed within first 
principles electronic structure calculationai3li3, in sev- 
eral cases exhibiting novel phenomena in nan nrn agnetism 
such as the existence of metamagnetic statesolij or inter- 
mixing effects between adatoms and the host surfacetLS. 
An accurate calculation of the total energy in terms of 
full potential or full charge density schemes made PQ^sir 
ble the investigations of the energetics of adatomstZTO. 
As compared to TB methods an obvious drawback of the 
embedded KKR technique is that, with respect to com- 
putational limitations, the number of the atoms in the 
cluster is restricted to about less than 100. Furthermore, 
the inclusion of structural relaxations is exceedingly dif- 
ficult. In order to circumvent these problems, a quasi-ab 
initio molecular-dynamics method has been employed by 
parameterizing interatomic potentials to the first jsrinci- 
ples KKR Green's function electronic structuresEj. On 
the level of a fully relativistic spin-polarized electron the- 
ory, recently, strongly enhanced orbital magnetism and 
MAE of adatoms and smaU-clusters on Ag and Au(100) 
surfaces have been reportedElrEj. 

From the mid-nineties on we carried out systematic 
investigations of the magnetism, in particular of the 
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MAE, of transition metal multilayer systems by using the 
fully relativistic spin-polarped-screened Korringa-Kohn- 
Rostoker (SKKR) methodE3t3. Specifically, within the 
single-site approximation, we explored the oscillatory be- 
havku: of the MAE of an Fe impurity buried in a Au(100) 
hostel The purpose of the present work is to extend 
these studies by including self-consistent effects (elec- 
tronic relaxations) of the host atoms in order to perform 
realistic investigations for magnetic clusters on metallic 
surfaces. For this very reason we make use of a real-space 
embedding technique in order to calculate the electronic 
structure of the cluster, and also to be able to treat the 
Poisson equation with appropriate boundary condition. 
Theoretical and computational details are given in Sec- 
tions H and III , respectively. In Section IV our results of 
the magnetic moments and the MAE of small planar Fe, 
Co and Ni clusters on Ag(100), as well as of the magnetic 
correlation between Fe or Co adatoms are presented. Fi- 
nally, in Section ^ we summarize and draw conclusions. 



II. THEORETICAL APPROACH 

Within multiple scattering theory the scattering path 
operator (SPO) matrix, r(E) = {z nm (E)} = {t™,{E)}, 
with Q and Q' being angular momentum indices and E 
being the energy, defined as 



r(E) = [t-\E)-G{E)Y 



(1) 



describes the full hierarchy of scattering effects between 
any two particular sites, n and m. In Eq. (|l|), t(E) = 
{t n {E) 6 nm } = {Vqq, (E) 6 nm } and G(E) = {G nm (E)} = 
{Gqqi (E)} denote the single-site i-matrices and the real- 
space structure constants, respectively. For more details, 
especially, how to calculate t r QQ,(E) within a fully rela- 
tivistic spin-polarized scheme, see, e.g., Ref.0. 

Assuming that a finite set of impurities interacts 
within a given finite range, we can select an environ- 
ment of impurities, C, containing also perturbed host 
atoms, such that outside C the potentials can be con- 
sidered to be identical with those of the unperturbed 
host. A particular cluster C can then be treated as 
perturbation of the host. In practice, we first calculate 
the SPO of the 2D translational invariant layered host, 
Th(k\\,E) = {t£ 9 (k|| , E)}, within the framework of the 
SKKR methodE3, where p and q denote layers and the 
k|| are vectors in the surface Brillouin zone (SBZ). The 
real-space SPO is then given by 

zr(E) = 7^— I e-^-^rlo^E^ , 



SBZ 



(2) 



where the atomic position vectors are decomposed as 
R m = Tj+Cp and R rl = Tj+c q with and Tj being 2D 
lattice vectors, c p and c q the so-called layer-generating 
vectors, and flssz is the unit area of the surface Brillouin 
zone. 



By replacing the i-matrices of the unperturbed host, 
th(E), with those of the cluster-atoms, tc(E), leads to 
the following Dyson like equation, 



t c (E) = r h (E) [I - {t- h \E) - t- c \E))T h {E)] ' 



(3) 



where tc{E) is the SPO-matrix corresponding to all sites 
in cluster C, from which in turn all corresponding local 
quantities, i.e., charge and magnetization densities, spin- 
and orbital moments, as well as the total energy can be 
calculated. Note, that Eq. (||) takes into account all scat- 
tering events both inside and outside the cluster. 

As usual, in order to perform self-consistent calcu- 
lations, the Poisson equation has to be solved repeat- 
edly. Partitioning the configurational space into non- 
overlapping cells centered at the atomic sites only the 
intracell part needs special care. By using a multipole 
expansion of the Coulomb potential, one arrives at the 
so-called Madelung potentials, which for a given site n 
of C will be denoted by V^f . Due to the additivity of the 
Poisson equation, V^ 1 can be decomposed into contribu- 
tions from atoms inside and outside of the cluster, V^ n 
and V^, respectively, 



v n M = v c M n 



V, 



M 



(4) 



where S denotes the ensemble of scatterers outside C. 
Clearly enough, assuming that the ensemble of sites in S 
is independent of C, implies that Vg 4 n has to be indepen- 
dent of the type of the atoms in C. Replacing the atoms 
in C by those of the unperturbed host, one can write 



t/M _ t/M,/i \ r M,h 
V S,n — V n ~ V C,n ' 

where V^ 1,h denotes the Madelung potential of sit 
the unperturbed, 2D translational invariant host 
V^ h stands for the Madelung potential of site n gener- 
ated by cluster C occupied by unperturbed host atoms. 
Substituting Eq. (|) into (|), yields 
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-n for 
3 and 



V, 



M 
C.n 



M,h 
C,n 



V, 



(6) 



which can be regarded as an embedding equation for the 
intercell potential resulting from the boundary condition 
set by the size of C. Obviously, by using Eq. (^) the 
problem of summing up the contributions to the intercell 
potential from region S is properly solved. 



III. COMPUTATIONAL DETAILS 

Self-consistent, fully relativistic calculations for se- 
lected, planar Fe, Co and Ni clusters on Ag(100) have 
been performed in the framework of the local spin- 
dcjisity approximation as parameterized by Vosko et 
alc3. In each case three different orientations for the 
magnetization were considered: along the z axis (normal 
to planes), as well as along the x and y axes (nearest 
neighbor directions in an fcc(lOO) plane). The potentials 
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were treated within the atomic sphere approximation 
(ASA). For the calculation of the i-matrices and for the 
multipole expansion of the charge densities, necessary to 
evaluate the Madelung potentials, a cut-off of £ max — 2 
was used. In order to perform the energy integrations, 16 
points on a semicircular contour in the complex energy 
plane were sampled according to an asymmetric Gaussian 
quadrature. Both, for the self-consistent calculation of 
the Ag(100) surface and for the evaluation of Eq. (0) we 
used 45 /c||-points in the irreducible wedge of the SBZ. 
For some restricted cases we checked the convergence of 
the results by increasing the number of fey-points to 210. 

In the present study we made no attempts to include 
lattice relaxation effects, therefore, the host and the clu&j 
ter sites refer to positions of an ideal fee parent latticed 
with the experimental Ag lattice constant (4.12 A). Three 
layers of self-consistently treated empty sites were used 
to represent the vacuum regions; the magnetic adatoms 
occupy sites in the first vacuum layer. As shown in Fig. |l|, 
we considered dimers and linear trimers oriented along 
the x a^4S, square-like tetramers, centered pentamers (as 
in Ref.E3), as well as a cluster arranged on positions of a 
3x3 square denoted in the following simply as 3 x 3 clus- 
ter. In Fig. |l|, for each particular cluster the equivalent 
atoms with respect to an orientation of the magnetization 
along the x or y axis are labelled by the same number. 
Note that for a magnetization aligned in the z direction, 
the atoms labelled by 2 and 3 in the pentamer and the 
3x3 cluster become equivalent. Up to a total of 67 sites, 
the clusters consisted of adatoms, several substrate Ag 
atoms and empty sites from neighboring shells. A sta- 
bility test of the local electronic and magnetic properties 
for a single Fe adatom with respect to the number of 
self-consistently treated neighboring shells is shown in 
Fig. 0. Although the calculated orbital moment of the 
Fe adatom shows a somewhat slower convergence than 
the valence charge and the spin moment, it is remark- 
able that considering only a first shell of neighbors this 
already yields values which differ by less than 1 % from 
the fully converged ones. 

IV. RESULTS AND DISCUSSION 
A. Spin and orbital moments 

Calculations for different orientations of the magneti- 
zation revealed that the spin moments are fairly insensi- 
tive to the direction of the magnetization, while for the 
orbital moments remarkably large anisotropy effects ap- 
ply, a phenomenon that will be discussed in the next Sec- 
tion. For a magnetization along the z axis, the calculated 
values of the spin and orbital moments for an adatom and 
selected clusters of Fe, Co and Ni on Ag(100) are listed 
in Table In there the position indices in a particu- 
lar cluster refer to the corresponding numbers in Fig. [j] 
and the number of nearest neighbors of magnetic atoms 
(coordination number, n c ) is also given. 

As compared to the corresponding monolayer values 
(3.15 hb for Fe and 2.03 hb for Co), the spin moment of a 



single adatom of Fe (3.39 hb) and Co (2.10 /x_b) is slightly 
increased. In the case of Fe clusters, the spin moments 
decrease monotonously with increasing n c . A slight de- 
viation from that behavior can be seen for the 3x3 clus- 
ter, where the atoms with n c —2 and 3 exhibit the same 
spin moment. For the central atom of the pentamer and, 
in particular, of the 3x3 cluster, the monolayer value is 
practically achieved. The above results compare fairly 
well to those of Cabria et alr3 and reflect a very short 
ranged magnetic correlation between the Fe atoms. 

The general tendency of decreasing spin moments with 
increasing n c is obvious also for the Co clusters up to 
the pentamer case. For the 3x3 cluster, however, just 
the opposite trend applies. Establishing a correlation 
between S z and n c for Co seems to be more ambiguous 
than for Fe, because the changes of the spin moment are 
much smaller in this case. Nevertheless, it is tempting to 
say that in the formation of the magnetic moment of Co, 
further off neighbors play a more significant role than in 
the case of Fe. 

In the case of an adatom and dimer of Ni we found 
no stable magnetic state. Quite contradictory, Cabria 
et alS3 reported a spin moment of about 0.5 [Lb for a 
Ni adatom on Ag(100). As the computational method 
of these authors is very similar to ours, it is at present 
puzzling what causes this discrepancy between the two 
calculations. Our result is, however, in line with the 
experiments of Beckmann and Bergmann who .found no 
magnetic moment for Ni adatoms on Au surfaceEll, which 
as a substrate is rathetsimilar to Ag. It should be noted, 
however, that in Ref.Ei the actual surface orientation is 
not specified. 

For clusters of Ni one can observe an opposite ten- 
dency as for Fe and Co: the spin moment enhances with 
increasing number of neighbors. This clearly can be seen 
from Table |j. Having in mind the calculated monolayer 
value (0.71 /is), our small cluster calculations indicate a 
fairly slow evolution of the spin moment of Ni with in- 
creasing cluster size, implying that the magnetism of Ni 
is subject to correlation effects on a much longer scale 
than in Fe or Co. 

Apparently, the orbital moments show a different, in 
fact, more complex behavior as the spin moments. For 
single adatoms of Fe and Co we found orbital moments 
enhanced by a factor of ~ 6 and ~ 4.5, respectively 
as compared to the monolayer values (0.14 /.is for Fe 
and 0.27 hb for Co). This is a direct consequence of 
the reduced crystal field splitting, being relatively large 
in monolayers, and, in particular, in corresponding bulk 
systemsEj. In spite of a qualitative agreement, our L z 
values for the adatoms are considerably larger than those 
calculated by Cabria et al!E3 (0.55/is for Fe and 0.76/is 
for Co). It should be noted, however, that by including 
orbital polarization effects (Himd's second rule)jJi terms 
of Brooks' parameterizationHEia, Nonas et alE3 found 
orbital moments for Fe and Co adatoms on Ag(100) close 
to the atomic limit (2.20/is for Fe and 2.57/xb for Co). 

For dimers of Fe and Co, the value of L z drops to 
about 40 % in magnitude as compared to a single adatom. 
The evolution of the orbital moment seems, however, to 
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decrease explicitly only for the—central atom of larger 
clusters. In a previous papera we showed that the 
(local) symmetry can be correlated with the magnetic 
anisotropy, i.e., with the quenching effect of the crystal 
field experienced by an atom. A single adatom and the 
central atom of the linear trimers, pentamers and the 
3x3 clusters exhibit well-defined rotational symmetry 
namely, C\, C2, C4 and C4, respectively. The corre- 
sponding values of L z , namely, 0.88 (13, 0.25 (13, 0.15 
fiB: and 0.12 [Ib for Fe, and 1.19 Hb, 0.49 [Xb, 0.25 \j,b 
and 0.23 [is for Co, nicely reflect the increasing rota- 
tional symmetry of the respective atoms. Although the 
outer magnetic atoms exhibit systematically larger or- 
bital moments than the central ones, even a qualitative 
correlation with the local environment (n c ) can hardly be 
stated. The orbital moment for the trimer of Ni is already 
close to the monolayer value (0.19 hb) but shows rather 
big fluctuations with respect to the size of the cluster and 
also to the positions of the individual atoms. 



B. Anisotropies of orbital moments and magnetic 
anisotropy energies 

By using the self-consistent potentials for a given ori- 
entation of the magnetization (along z), we calculated 



magnetic anisotrepy energies by means of the magnetic 
force theorerr£3'E£l as differences of band-energies, 



E\ YT — Eh-, and AE„-, — Eh-,, — Eh-, 



(7) 



For a particular orientation a, the band-energy is ob- 
tained as a sum of contributions from all atoms in the 
cluster. 



E b . a = 'J2El. a (a = x,y,z), 



iec 



El 



de (e - e F ) <(e) , 



(8) 



(9) 



where e F is the Fermi energy of the substrate, es is the 
bottom of the valence band and n z a (e) is the density of 
states for atom i. Clearly, the above formalism allows us 
to define the MAE as a sum of atom-like contributions, 
which facilitates to trace its spatial distribution in the 
cluster. 

The anisotropies of the orbital moments and the con- 
tributions of the individual magnetic atoms to the MAE 
are displayed in Tables ||, |n| and IV for Fe, Co and Ni 
clusters, respectively. In addition, the total MAE per 
magnetic atoms of the clusters including the neighbor- 
hood is also given. Although the dominating contribu- 
tions to the MAE arise from the magnetic species, the 
environment, in particular, the Ag atoms and the empty 
sites within the first shell add a remarkable amount to the 
MAE. However, due to the weak polarization of the Ag 
atoms, we obtained a fast convergence of the total MAE 
with respect to the size of the cluster (environment). 



As can be inferred from the corresponding positive val- 
ues of the MAE in Tables || and III, single adatoms of 
Fe and Co exhibit a magnetization oriented perpendic- 
ular to the surface. This again is in perfect agreement 
with the experiments of Beckmann and Bergmanr£3. As 
compared with the monolayer case (0.47 meV), the MAE 
of an Fe adatom (5.61 meV) is enhanced by— a factor of 
twelve. Contrary to our results, Cabria et alxin predicted 
in-plane magnetism (AE X - Z = -0.98 meV) for an Fe 
adatom on Ag(100), and perpendicular magnetism for 
Co, albeit with a much larger anisotropy energy (> 7 
meV) than ours (4.36 meV). It should be stressed at 
this point, that our calculations are consistent with a 
qualitative rule, valid for transition metals with a more 
than half-filled ci-band and based on simple, pecturba- 
tive phenomenological or tight -binding reasoningci the 
direction, along which the orbital moment is the largest, 
is energetically favored. 

As can be seen from Table ||, perpendicular mag- 
netism is characteristic for all Fe clusters considered. For 
the dimer and the trimer we observe a small in-plane 
anisotropy with preference of the x axis, i.e., in the di- 
rection of the Fe-Fe bonds. In agreement with the reduc- 
tion of the orbital moment, as discussed in the previous 
section, the contribution of the central atom to the MAE 
for the trimer, the pentamer and the 3x3 cluster rapidly 
decreases, being even less than the monolayer value in 
the case of the 3x3 cluster. The outer atoms in the pen- 
tamer and in the 3x3 cluster can add considerably more 
to the MAE than the central atom. As a consequence, 
the average MAE strongly fluctuates with increasing size 
of the magnetic cluster and shows a very slow tendency 
to converge to the MAE of an Fe monolayer on Ag(100). 
Such a complicated behavior of the MAE with respect t< 
the cluster size has also been found by Guirado-Lopezl 
for free-standing fee transition metal clusters. 

In comparison to an adatom, for a Co dimer AE X - Z 
drops to a large negative value (-3.49 meV/per Co), while 
AEy- z remains slightly positive (0.76 meV/per Co), im- 
plying that a Co dimer favors the x (in-plane) direction 
of the magnetization and also experiences a strong in- 
plane anisotropy. The strong tendency of Co clusters to 
in-plane magnetization pertains to larger clusters and is 
characteristic also for a Co monolayer (AE X - Z = -1.31 
meV). The atom-like resolution of the MAE indicates, 
that this tendency is driven by nearest-neighbor Co-Co 
interactions. An explanation of this effect in terms of 
perturbation theory and symmetry resolved densities of 
states can be found in Refs.ElTta. As an unexpected 
consequence, the contribution to the MAE of the cen- 
tral atom in the cluster can be larger than that of some 
outer atoms. Quite obviously, the MAE of the central 
atom of the trimer, the pentamer and the 3x3 cluster, 
-9.06 meV, -2.46 meV, and -1.86 meV, respectively, fall 
monotonously off to the monolayer value, whereas the av- 
erage MAE possesses a much more complicated evolution 
also in this case. 

With exception of the tetramer, for which we found 
a MAE close to zero, all Ni clusters prefer an in-plane 
magnetization. The in-plane anisotropy, seen from Ta- 
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ble IV for the trimer, but also from the atom-like contri- 



butions for the larger clusters, is, however, smaller than 
in the case of Co. Again the complicated nature of the 
magnetism of Ni shows up, in particular, for the 3x3 
cluster: while the contribution of the central atom to the 
MAE almost vanishes, those of the outer atoms oscillate 
in magnitude. Considering the MAE of a Ni monolayer 
on Ag(100) (-2.23 meV), no straightforward connection 
with the magnetic anisotropy properties of small clusters 
can be traced. 



to a parallel configuration, the corresponding values of 
S z and L z in Table | can be questioned. Therefore, for 
this particular case we performed self-consistent calcula- 
tions also for the antiparallel alignment. Assuringly, for 
both Fe and Co, we obtained the same value of S z and 
L z within 1% relative accuracy as in the case of a parallel 
alignment. For larger distances we observe ferromagnetic 
coupling, which virtually vanishes for d > 5a, implying 
a very weak, short ranged exchange interaction between 
the adatoms of Fe and Co induced by the Ag host. 



C. Magnetic interaction between adatoms 

Interactions between magnetic nanoclusters are of 
great importance for technological applications. Clearly 
enough the most important questions are (i) what is the 
magnetic structure of the individual entities, (ii) of what 
nature (strength, range, etc.) is the coupling between 
them, and (Hi) what influences the magnetic orientation 
of these entities relative to each other. In this section we 
present a preliminary study in this field by investigating 
the interaction of two Fe or Co adatoms on Ag(100). 

We first performed self-consistent calculations for two 
adatoms by varying the distance d between them from a 
to 5a along the x direction, where a is the 2D lattice con- 
stant and keeping the orientation of the magnetizations 
parallel to each other (along the z axis). The calculated 
spin and orbital moments of the (coupled) adatoms are 
shown in Fig. |[ Note that the distance a refers to the 
bondlength in dimers. As can be seen from Fig. both 
for Fe and Co the values of S z and L z rapidly converge 
to the respective single adatom value. 

Next we calculated the exchange-coupling energy, 
AEx , between the two adatoms by taking the energy dif- 
ference between a parallel (TT) and an antiparallel (||) 
orientation of the two adatoms, 

AE X = E 6 (TT) - E b (U) ■ (10) 

using, however, the self-consistent potentials for the par- 
allel configuration. We are aware of the fact that, due to 
the lack of self-consistency in the antiparallel case, for 
near adatoms this approach might be quite poor. We be- 
lieve, however, that this approximation provides a good 
estimate of the sign and the magnitude of the interaction. 

The calculated AEx is shown in Fig. ^ for Fe and Co 
as a function of the distance d between the two adatoms. 
Apparently, for d = a in both cases a strong, ferromag- 
netic nearest-neighbor exchange-coupling between these 
two atoms applies, with an interaction energy somewhat 
larger for Fe than for Co. As the two adatoms are adja- 
cent in this case, this strong coupling can be attributed 
to a direct exchange mechanism. Increasing the separa- 
tion between the two adatoms, AEx rapidly decreases. 
For d = 2a it changes sign, i.e., the coupling becomes 
antiferromagnetic. Since for an antiparallel alignment of 
the spin moments of the two adatoms, lying close to each 
other, the electronic structure and the magnetic moments 
might be expected to differ to some extent as compared 



V. SUMMARY AND CONCLUSIONS 

By using a real-space embedding technique based on 
the Korringa-Kohn-Rostoker Green's function method, 
we have performed fully relativistic, self-consistent cal- 
culations for adatoms and small clusters of Fe, Co and 
Ni on Ag(100). Due to the decreased coordination of the 
magnetic atoms, we obtained slightly enhanced spin mo- 
ments for adatoms and small clusters of Fe and Co and 
found that the spin moments are already close to the 
monolayer values_£pr a cluster of 9 atoms. In agreement 
with experimentsEB the adatoms and dimers of Ni turned 
out to be nonmagnetic, while the spin moments in larger 
Ni clusters indicated a complex formation of magnetism. 
In connection with strongly enhanced orbital moments, 
for Fe and Co adatoms we revealed an unusually strong 
tendency to perpendicular magnetism. The perpendic- 
ular magnetism persisted also for Fe clusters of increas- 
ing size, whereby the atom-like contributions showed an 
oscillating behavior depending mainly on the local rota- 
tional symmetry. The preferred orientation for clusters of 
Co and Ni obtained was in-plane. In addition, we inves- 
tigated the magnetic coupling between two adatoms of Fe 
or Co, for which we established a 'good local-moment' 
behavior. In terms of calculated exchange-coupling en- 
ergies, the dimers show a strong ferromagnetic coupling, 
which immediately drops two orders of magnitude with 
increasing distance between the two adatoms, indicating 
a weak, indirect coupling between them. 

The main outcome of the present paper is that by per- 
forming first principles calculations, not only the qual- 
itative trends of small cluster magnetism of transition 
metals, but even quantitative results can be obtained 
which in turn can directly be compared with experiments. 
By using a parallelized version of our computer code the 
number of atoms treated in the cluster can be easily in- 
creased to some hundreds. This is encouraging to extend 
the present calculations to larger nanostructures (mag- 
netic wires, dots, corrals etc.) currently being in the 
very focus of technological applications. 
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TABLE I. Calculated spin moments (S z ) and orbital mo- 
ments (I/ z ), in units of \ib, for small clusters of Fe, Co and Ni 
on Ag(100) with magnetization perpendicular to the surface 
(z). For each position in a particular cluster (see Fig. |l|), n c 
refers to the number of the neighboring magnetic (Fe, Co, Ni) 
atoms. 



Fc 



Co 



Ni 



Cluster 


position 


n c 


s z 




s z 








adatom 







3.39 


0.88 


2.10 


1.19 






dimer 


1 


1 


3.31 


0.32 


2.09 


0.49 






trimer 


1 


2 


3.29 


0.25 


2.07 


0.45 


0.77 


0.21 




2 


1 


3.33 


0.44 


2.06 


0.49 


0.70 


0.23 


tetramer 


1 


2 


3.26 


0.18 


2.08 


0.32 


0.76 


0.28 


pentamer 


1 


4 


3.13 


0.15 


2.01 


0.25 


0.76 


0.12 




2 


1 


3.35 


0.37 


2.10 


0.59 


0.71 


0.33 




3 


1 


3.35 


0.37 


2.10 


0.59 


0.71 


0.33 


3x3 cluster 


1 


4 


3.15 


0.12 


2.06 


0.23 


0.79 


0.24 




2 


3 


3.23 


0.16 


2.04 


0.30 


0.71 


0.20 




3 


3 


3.23 


0.16 


2.04 


0.30 


0.71 


0.20 




4 


2 


3.23 


0.33 


2.00 


0.29 


0.63 


0.19 



G 



TABLE II. Calculated orbital moment anisotropies (AL), 
in units of /is , and contributions of the Fe atoms to the MAE, 
AE, in units of meV, for small clusters of Fe on Ag(100). 
For each cluster, the total MAE per Fe atom of the cluster 
including the neighborhood is also given in parentheses. 



Cluster 

v_yl LAO IC1 






AT, 


AE 


AT, 


AE 


adatom 







-0.37 


5.07 


-0.37 


5.07 










(5.61) 




(5.61) 


dimer 


1 


1 


-0.12 


2.14 


-0.11 


1.66 










(2.30) 




(1.83) 


trimer 


1 


2 


-0.12 


1.93 


-0.08 


0.93 




2 


1 


-0.16 


2.83 


-0.15 


2.39 










(2.72) 




(2.13) 


tetramer 


1 


2 


-0.02 


0.50 


-0.02 


0.50 










(0.54) 




(0.54) 


pentamer 


1 


4 


-0.03 


0.49 


-0.03 


0.49 




2 


1 


-0.03 


0.92 


-0.08 


0.85 




3 


1 


-0.08 


0.85 


-0.03 


0.92 










(0.90) 




(0.90) 


3x3 cluster 


1 


4 


0.00 


0.23 


0.00 


0.23 




2 


3 


-0.02 


0.43 


-0.01 


0.84 




3 


3 


-0.01 


0.84 


-0.02 


0.43 




4 


2 


-0.13 


1.86 


-0.13 


1.86 










(1.20) 




(1.20) 


TABLE III. 


As in 


r 

Table |l 


l| for Co clusters. 




Cluster 


position 


Tic 


A J 


AEx-z 


ALy — Z 


A F 


adatom 







-0.26 


4.20 


-0.26 


4.20 










(4.36) 




(4.36) 


dimer 


1 


1 


0.15 


-3.50 


-0.01 


0.67 










(-3.49) 




(0.76) 


trimer 


1 


2 


0.40 


-9.06 


-0.02 


-0.11 




2 


1 


0.34 


-6.29 


0.05 


-0.04 










(-7.44) 




(-0.01) 


tetramer 


1 


2 


0.15 


-2.29 


0.15 


-2.29 










(-2.37) 




(-2.37) 


pentamer 


1 


4 


0.12 


-2.46 


0.12 


-2.46 




2 


1 


0.21 


-4.16 


-0.01 


-0.03 




3 


1 


-0.01 


-0.03 


0.21 


-4.16 










(-2.22) 




(-2.22) 


3x3 cluster 


1 


4 


0.13 


-1.86 


0.13 


-1.86 




2 


3 


0.10 


-1.56 


0.18 


-2.96 




3 


3 


0.18 


-2.96 


0.10 


-1.56 




4 


2 


0.16 


-2.60 


0.16 


-2.60 










(-2.45) 




(-2.45) 



TABLE IV. 


As 


in Table [I 


J] for Ni 


clusters. 




Cluster 


position 


Tl r 


AL . 


L -* J — / x — z 


A L„ . 


l—±Uy — z 


trimer 


1 


2 


0.19 


-6.12 


0.18 


-1.38 




2 


1 


0.11 


-3.72 


0.08 


-1.00 














(-1.13) 


tetramer 


1 


o 

I 


-O.UO 




-O.UO 


U.lH 










(0.10) 




(0.10) 


pentamer 


1 


A 

4 


ft 1 r 
U.10 


-z.zo 


u.io 


-Z.Zo 




2 


1 


a a i 


-1.04 


U.UO 






Q 
O 


1 


one; 
U.UO 


-U.Dy 


A A1 
U.U1 


-1.D4 










(-1.41) 




1 1 A~\\ 

(-1.41) 


3x3 cluster 


1 


4 


-0.06 


-0.02 


-0.06 


-0.02 




2 


3 


0.02 


-0.75 


0.06 


-2.00 




3 


3 


0.06 


-2.00 


0.02 


-0.75 




4 


2 


0.05 


-1.21 


0.05 


-1.21 










(-1.17) 




(-1.17) 



7 



®© 

0© ©©© ®® 

dimer trimer 



tetramer 



© 
@®@ 
© 

pentamer 



®@® 
®Q® 
®@® 

3x3 cluster 



FIG. 1. Sketch of the planar clusters considered. For an 
orientation of the magnetization along the x or y axis, the 
equivalent atoms in a cluster are labelled by the same number. 
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FIG. 2. Calculated number of valence electrons (N va i), spin 
moment (S z ) and orbital moment (L z ) of a single Fe adatom 
on a Ag(100) surface as a function of the number of the 
self-consistently treated atomic shells around the Fe atom. 
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FIG. 3. Calculated spin and orbital moments of two 
adatoms of Fe or Co on Ag(100) as a function of their dis- 
tance d measured in units of the 2D lattice constant a. 
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FIG. 4. Calculated exchange coupling energy, AEx, be- 
tween two adatoms of Fe or Co on Ag(100) as a function of 
the distance d measured in units of the 2D lattice constant a. 
The insets show the range 2a < d < 5a on a blown up scale. 
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